Angular dependent planar metamagnetism in 
the hexagonal compounds TbPtIn and TmAgGe 
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Detailed magnetization measurements, M(T,H,S), were performed on single crystals of TbPtIn 
and TmAgGe (both members of the hexagonal Fe2P/ZrNiAl structure type), for the magnetic field 
If applied perpendicular to the crystallographic c axis. These data allowed us to identify, for each 
compound, the easy-axes for the magnetization, which coincided with high symmetry directions 
([120] for TbPtIn and [110] for TmAgGe). For fixed orientations of the field along each of the 
two six-fold symmetry axes, a number of magnetically ordered phases is being revealed by M(H,T) 
measurements below Tjv. Moreover, T ~ 2 K, M(H)|e measurements for both compounds (with H 
applied parallel to the basal plane), as well as T — 20 K data for TbPtIn, reveal five metamagnetic 
transitions with simple angular dependencies: Hci,j ~ 1/ cos(6 ± ip), where (p = O" or 60". The 
high field magnetization state varies with 9 like 2/3 * jj,sat(R^'^) * cos 9, and corresponds to a crystal 
field limited saturated paramagnetic, CL-SPM, state. Analysis of these data allowed us to model 
the angular dependence of the locaUy saturated magnetizations Msat and critical fields He with a 
three coplanar Ising-like model, in which the magnetic moments are assumed to be parallel to three 
adjacent easy axes. Furthermore, net distributions of moments were inferred based on the measured 
data and the proposed model. 

PACS numbers: 75.25.-(-z; 75.10.-b; 75.30.Gw; 75.30.Kz; 75.50.Ee 
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I. INTRODUCTION 



Numerous studies of angular dependent local moment 
metamagnetism show that simple geometric relationships 
exist between the critical fields of the metamagnetic 
phase transitions, and the angle that the applied field 
makes with the corresponding easy axis. One simple case 
is that of TbNi2Ge2^, a tetragonal compound with Tb 
ions in tetragonal point symmetry, where, at low tem- 
peratures, the crystalline electric field (CEF) anisotropy 
confines the local moments along the c ([001]) crystal- 
lographic axis (Ising-like system). Consequently, several 
metamagnetic transitions are observed, with the critical 
field values He having a 1/ cos 9 dependence on the an- 
gle between the applied field H and the c-axis. A more 
complex situation is encountered when the magnetic mo- 
ments are allowed more degrees of freedom , i.e. when the 
CEF anisotropy constrains them to an easy plane. This 
is the case in the tetragonal compounds RNi2B2C— iii, 
with R = Tb - Er, and RAgSb2^ for R = Dy. The R 
ions are again in tetragonal point symmetry and the lo- 
cal moments are confined to four equivalent [110] or [100] 
crystallographic directions; thus the angular dependent 
magnetization measurements, when the field is applied 
in the basal plane, reveal the four-fold anisotropy of the 
longitudinal magnetization that reflects the symmetry of 
the unit cell. The angular dependencies of the locally 
saturated magnetizations Msat and the critical fields He 
could be treated by simple analysis, and plausible net 
distribution of moments could be inferred for each meta- 
magnetic phasc"^'^. Kalatsky and Pokrovsky^ elaborated 
the four-position clock model, which agrees well with the 
observed metamagnetism in HoNi2B2C and DyAgSb2. 



Initially thought to be improbable, angular dependent 
metamagnetism in extremely planar systems is now an 
accepted and understandable event in tetragonal com- 
pounds; but this has not yet been well studied in hexago- 
nal compounds, and this is the motivation for the present 
analysis. Recently we reported the presence of meta- 
magnetism in a hexagonal class of materials, i.e. the 
RAgGe compounds, for R = Tb - Tm^. They crystal- 
lize in the ZrNiAl-typc structure, an ordered variant of 
the hexagonal Fe2P family. In this structure, there is a 
unique rare earth site in the unit cell, with the rare earth 
ions occupying equivalent 3g positions with orthorhom- 
bic point symmetry (fig.l). The anisotropic susceptibility 
was found to be axial for TbAgGe and progressed towards 
extremely planar for TmAgGe. The physical properties 
of the TmAgGe compound (antiferromagnetic ordering 
in the ground state, extremely planar anisotropy, meta- 
magnetism when the applied field was perpendicular to 
the c-axis) are very similar to those of the isostructural 
TbPtIn system, previously only known in polycrystalline 
form®'^, making both compounds good candidates for a 
study of the angular dependent metamagnetism in hexag- 
onal systems. In both TmAgGe and TbPtIn the R ions 
occupy the same site, whereas both ligands are differ- 
ent. Having two such systems will allow us to show 
that the behavior we find is not specific to one com- 
pound, but a more general result associated with this 
structure, or perhaps with the orthorhombic point sym- 
metry in a hexagonal unit cell. As we shall see, the Neel 
temperature is much higher for TbPtIn (46.0 K) than for 
TmAgGe (4.2 K), with the former also showing a possible 
spin reorientation transition at a fairly high temperature 
(T,„ = 27.4 K). This latter transition in TbPtIn was 
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missed by the measurements on polycrystalline samples 
where even the nature of the magnetically ordered state 
below ~ 50 K was not identified. Thus we can once again 
emphasize the advantage of analysis on single crystals 
rather than on polycrystalline samples. 

This paper is organized as follows: after a brief descrip- 
tion of the experimental methods, we present the M(T,H) 
data on TbPtIn, emphasizing the complexity of its or- 
dered state, for field orientations along the two in-plane, 
high symmetry directions; this is followed by the M(H,6') 
measurements at low temperature, from which the values 
of the locally saturated magnetizations and critical fields 
as a function of 9 were extracted. Next, we will introduce 
a model for the net distribution of moments, which we 
subsequently use to calculate the expected, locally satu- 
rated magnetizations Mj and the critical fields (for the 
transition from state i to state j) Hdj as functions of 
0. This will be followed by a comparison of how the cal- 
culated and experimental Mj and Hd^j values vary with 
the angle between the applied field H and the easy axis 
for this compound. 

Similar measurements performed on TmAgGe will 
then be shown, leading to the corresponding experimen- 
tal Mj{9) and Hcij{0) phase diagrams, which we will 
again map onto the appropriate model calculation. 

We will also analyze similar M{H, 9) data of TbPtIn 
for T — 20 K, and use the same model in order to char- 
acterize a magnetic phase present only at higher temper- 
atures. 

Finally, we will summarize our key results and also 
indicate how our model can be generalized to a variety 
of possible point symmetries in tetragonal and hexagonal 
systems. 



II. EXPERIMENTAL METHODS 

Using an initial concentration of Tbo.05Pto.05Ino.90j 
clean hexagonal rods of TbPtIn were grown; this In-rich 
self-flux was used because of its low-melting tempera- 
ture and because it introduces no new elements into the 
melt. The constituent elements were placed in alumina 
crucibles and sealed in quartz tubes under partial argon 
pressure; subsequently they were heated up to 1200° C, 
and then cooled down to 800° C over 48 h. Removal of 
the excess liquid solution revealed hexagonal rods with 
the c-axis along the axis of the rod. The crystal struc- 
ture was confirmed by room temperature powder x-ray 
diffraction measurements, using Cu Kq, radiation, from 
which the lattice parameters a = (7.55 ± 0.01) A and c 
= (3.88 ± 0.01) A were also obtained. It is worth noting 
that the above temperature profile yielded in most cases 
well-formed, fairly good crystals (residual resistivity ratio 
RRR = 5.2); however, some of the crystals had hollow 
channels in the center, sometimes with flux inclosures, 
and the quality of the resulting hexagonal rods seems to 
only deteriorate with slower cooling rates {i.e. 400° C / 
lOOh). Single crystals of TmAgGe were also flux-grown 



out of Ag-Ge self flux, as described in detail in Morosan 
et al7 . 

Magnetic measurements were performed in a Quan- 
tum Design Magnetic Properties Measurement System 
(MPMS) SQUID magnetometer, with a specially mod- 
ified sample holder to allow the rotation of the sam- 
ple so that its c-axis stays perpendicular to the applied 
magnetic field. Additional measurements up to 140 kG 
were also taken, using an extraction magnetometer in a 
Quantum Design Physical Properties Measurement Sys- 
tem (PPMS). In order to avoid torque on the rotator due 
to the extreme anisotropy of the samples , small mass 
samples {i.e. m = (0.40 ± 0.05) mg in the case of 
TbPtIn, and m= (0.54 ± 0.05) mg for TmAgGe) were 
used for the angular dependent measurements. To cor- 
rect for the large weighing errors that would result from 
such small masses, we used a 13.55 mg TbPtIn sample, 
and a 2.40 mg TmAgGe sample respectively, to measure 
M(H) curves for 6* = 0° and 30° (with the angle 9 mea- 
sured between the applied field H and the easy axis, as 
described below). Whereas some errors are still intro- 
duced by this calibration due to the manual orientation 
of the large samples, the angular accuracy with which we 
are able to orient them is probably within 10°. The var- 
ious M(H, 9) measurements on the smaller samples were 
than calibrated using the data for the larger mass pieces; 
this is believed to be very accurate given that the angular 
uncertainty in the rotator is less than 1°. There is addi- 
tional uncertainty introduced by possible misalignment 
of the small piece with the c-axis exactly perpendicular 
to the applied field (due to the construction of the sample 
holder, this misalignment shouldn't be more than 10°). 

Laue experiments have been performed on well-formed 
single crystals. They confirmed that the as-grown 
facettes of the hexagonal rods are alternatingly paral- 
lel to the (101) and (Oil) crystallographic planes (the 
macroscopic facettes are aligned with the hexagonal unit 
cell axes [100] or [010]). 



III. TbPtIn 

Fig. 2 shows the inverse magnetic susceptibility for 
TbPtIn, as well as the low-temperature magnetization 
in the inset. This compound appears to order anti- 
ferromagnetically below T/v — 46.0 K, with what 
is probably a spin-reorientation or a commensurate-to- 
incommensurate transition around T,„ — 27.4 K, as 
indicated by the peaks in the d{M * T/H)/dT plot^o in 
Fig. 3a. In the determination of these temperature values, 
M/H represents the polycrystalline average susceptibil- 
ity Xave, calculated as 

1/ — 1 1/ _L 2 , , 

Xave — 3 AC "r 3 X.ab 

or, when measurements have been done along all three 
salient directions, 

Xave = 3 (X[001] + X[100] + X[120])- 

The above temperatures are further confirmed by the 
Cp{T) and p{T) (j||a6) data, also shown in Fig. 3b and c. 
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In the ordered state, as well as in the paramagnetic 
state up to ~ 100 K, the susceptibility is extremely 
anisotropic, with the local Tb moments confined to the 
a&-plane at low temperatures. To check the origin of 
the anisotropy in the paramagnetic state, single crys- 
tals of YPtIn were grown, with a small number of non- 
magnetic Y'^+ ions substituted with magnetic Tb'^^ ions. 
The extremely anisotropic susceptibility and magneti- 
zation of the diluted compound demonstrate that this 
is single-ion anisotropy associated with the CEF split- 
ting of the Hund's rule ground state J multiplet. From 
the Curie- Weiss effective moment, as determined from 
the average inverse susceptibility (Fig. 4), the concen- 
tration X of the diluted compound (Tb2:Yi_2;)PtIn is 
X — 0.024, whereas the high field magnetization equals 
~ 5.6 /XB/Tb'^+. (If we assume the high field magneti- 
zation to be ~ 6.25/iB/Tb'^^, as expected based on the 
model described below, the resulting concentration will 
be a; = 0.019). 

For temperatures higher than 150 K, Curie- Weiss be- 
havior of the pure TbPtIn compound can be inferred 
(Fig. 2) from the linear inverse susceptibilities, resulting 
in anisotropic Weiss temperatures Qab — 39.0 K and 
&c — 23.3 K. The polycrystalline average susceptibility 
Xave yielded an effective moment ^e// = 9.74/iB/Tb^+, 
very close to the theoretical value 9.72 /xb; the corre- 
sponding Weiss temperature is Qave = 33. 3K. 

The field dependent magnetization measurements 
shown in Fig. 5a not only confirm the in-plane/out-of- 
plane anisotropy observed in the ordered state, but 
also indicate anisotropic magnetization within the basal 
plane. Moreover, several metamagnetic transitions can 
be seen for the field parallel to the a6-plane, for fields 
up to 140 kG in the AI (H) data, and for the magne- 
toresistance measurements up to 90 kG in Fig. 5b. The 
geometry of the crystals led to more uncertainty in orient- 
ing the resistance pieces than the ones for magnetization 
measurements; therefore we can infer the approximate 
orientation of the magnetoresistance sample with respect 
to the applied field, by comparing the Ap{H) / p{0) data 
(Fig. 5b) with the M{H) curves for H\\ab (Fig. 5a): since 
various features in the magnetoresistance measurements 
occur closer to the critical fields in Af ([120]), we can as- 
sume that the field was almost parallel to the [120] di- 
rection. (The sharp drop in p{H) below ^ 1 kG (Fig. 5b) 
is very likely due to superconductivity of residual In flux 
on the surface of our resistance bar). For H along the 
c-axis, the magnetization increases almost linearly with 
increasing field (Fig. 5a), while staying far smaller than 

Mab- 

At the highest applied field {H — 140 kG), the mag- 
netization values for the three shown orientations are 
M([110]) = 5.86Mi3/Tb, Af([120]) = GAbfis/Th 
and M([001]) = 0.92^B/Tb. Whereas the extreme 
planar anisotropy of TbPtIn and the anisotropy within 
the a&-plane recommended this compound for a study 
of the angular dependent metamagnetism, the fact that 
the magnetization values were smaller in all three direc- 



tions than the calculated figat = ^fJ-B for T&^+ ions is 
somewhat intriguing. One plausible explanation for the 
low magnetization values would be the existence of more 
metamagnetic transitions for fields unaccessible with our 
measurement systems {i.e. above 140 kG). Another pos- 
sibility is that an additional energy scale (such as CEF 
splitting) exists, that confines the three local moments 
to three distinct, non-coUinear, in-plane orientations. As 
shall be shown below, we believe the latter to be the more 
likely scenario. 

In order to determine the easy axes of the system, 
we continuously rotated a small piece of the diluted 
sample (Tba;Yi_a;)PtIn {x w 0.02) in an applied field 
H = 55 kG (perpendicular to c); for constant tem- 
perature T = 2 K, the corresponding magnetization 
measurement, shown as open symbols in Fig. 6, roughly 
follows a cosO dependence (solid line) around the clos- 
est [120] -equivalent directions {i.e. 9 — 60° * n, 
n = integer), where the maxima occur. As we have 
seen before, the concentration of the diluted sample is 
X = (0.0215 ± 0.0025); given this uncertainty, the ab- 
solute value of the magnetization M{9) in Fig. 6 could not 
be determined, and thus we scaled the data to the max- 
imum value, M{0 = 0°). Similar behavior appears in 
the pure TbPtIn compound (full symbols in Fig. 6), where 
the measured data have also been scaled to their corre- 
sponding maximum value at = 0°. We can conclude 
that the easy axes of the TbPtIn system are the [120] 
directions. However dramatic departures from the cos 9 
angular dependence can be noticed. The magnetization 
for TbPtIn, indicative of strong interactions between the 
local moments, is also consistent with various metamag- 
netic states crossing the H — 55 kG line at different 
angles. Based on the above data, we will consider the 
easy axes to be the [120]-equivalent directions, and the 
angle 9 will be measured from the closest easy axis. 

In order to get an idea about the various metamagnetic 
states in this compound, we first explored changes of the 
critical fields and temperatures for two fixed orientations. 
The corresponding M{T) \H,e and M{H) [y.e measure- 
ments shown in Fig. 7-8 have been used to determine the 
H—T phase diagrams for the two in-plane high symmetry 
directions, -ff||[120] and i?||[110] respectively. As illus- 
trated in the insets in these figures, the points in these 
phase diagrams have been determined from local max- 
ima in d{M * T)/dT (fuU circles in fig.9) for fixed fields, 
and in dM/dH (open circles in fig.9) for various M{H) 
isotherms. Even though in Fisher et ali^ the maxima in 
d{M * T) /dT criterion is described only for antiferromag- 
netic systems, we apply it here not only for the AF state, 
but also for high-field states, where the magnetization has 
a net ferromagnetic component. We are confident that 
small errors are thus introduced, given the consistency of 
the critical field and temperature values obtained from 
both d{M*T)/dT and dM/dH derivatives (fuU and open 
symbols respectively in Fig.9). Given that the transition 
peaks were broad for some field and temperature values, 
we used Lorentzian fits of the corresponding derivatives 
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(thick lines in fig. 7 and 8, insets) to determine the critical 
values and T^. 

The resulting H — T phase diagrams for if||[120] 
(fig. 9a) and i7||[110] (fig. 9c) arc qualitatively similar, at 
low temperatures and low fields showing the metamag- 
netic states already seen in the M(H) data in Fig. 5. For 
i7||[120], fig. 9a shows that the antiferromagnetic groimd 
state persists up to about 20 kG, followed by a small 
intermediate state Mi (between ~ 20 kG and 28 kG) 
and a higher field state M2 up to ~ 54 kG; as field is 
being further increased, the paramagnetic PM state is 
reached, as already indicated in fig. 5a by the horizontal 
plateaux measured up to 140 kG. (At low temperatures, 
this is a crystal-field limited saturated paramagnetic CL- 
SPM state, in which, as discussed below, all moments are 
assumed to be in their 'up' positions, while still confined 
by the strong CEF energy to three distinct, non-collinear 
directions within the basal plane.) When moving up in 
temperature at low fields, we find the antiferromagnetic 
ground state to extend up to ^ 27 A K, whereas the mag- 
netic ordered state persists up to ~ 46.0 K; both transi- 
tion temperatures have been already observed (fig. 2 and 
3). The Ml phase exists below ^ 5.0 K, after which, for 
a limited temperature range (5.0 K < T < 15.0 K), 
there is a direct transition from the AF to the M2 state. 
Between 15.0 K and 27.4 K, or 18.4 kG and 2.0 kG respec- 
tively, another intermediate phase, M4, forms. The inset 
in fig. 8a represents an example of two isothermal cuts of 
the [120] phase diagram in fig. 9a: the phase boimdaries 
of the bubble-like phase Af4 are very close in field at con- 
stant T, and so the lower peak in the T = 20.0 K isother- 
mal derivative dM/dH is poorly defined; the two higher 
peaks are fairly sharp, similar to the one defining the M2 
to PM transition in the T ~ 30.0 K isotherm. However, 
as we move down in field along the latter isotherm, one 
broad peak around 15.4 K may indicate the crossing of 
another almost horizontal phase boundary (leading into 
the M5 state), and thus hard to identify in d{M*T)/dT. 
One more peak at He — 3.4 kG is indicative of possibly 
another phase AF' existing below this field, between 28 K 
and 46.0 K. This is consistent with an antiferromagnetic 
ordered state below 46.0 K, with an inconunensurate- 
commensurate transition around 28 K which frequently 
occurs in intermetallic compounds. 

For /fllillO] (Fig.9c) the iJ-T phase diagram is fairly 
similar, with only a few differences: a far less distinct 
Ml phase and a lower upper-boundary for the M2 re- 
gion. However, the most notable difference is a new high 
field phase, M3, whose upper boundary is determined 
by the points indicated with small arrows in fig. 7b. (As 
the field is being increased towards 140 kG, this line be- 
comes almost vertical, making it difficult to identify also 
in M{H)\t measurements.) 

In order to see how the H — T phase diagram 
evolves from if||[120] to if||[110], we collected compa- 
rable M{T)\h and M{H)\t data for an intermediate ori- 
entation of the applied field (approximately 12*^ from the 
easy axis [120]). Fig. 10a shows these M{T) curves for 



H = 1 — 70 kG, with the small arrows indicating the 
highest-T transition at each field value, as determined 
from the d{M * T)/dT maxima. Thus, the upper-most 
phase boundary in Fig. 9b, representing the phase dia- 
gram for this intermediate position, can be followed in 
field np to H = 65 kG (full symbols). For this orien- 
tation, we can also identify this line in the M{H) data, 
and an example is shown in Fig. 10b for the T = 10 K 
M{H) isotherm and its dM/dH derivative. Overall, the 
features common to both Fig. 9a and c are also present 
in Fig. 9b; moreover, in going from the [110] to the [120] 
direction, the M3 phase is being compressed, such that in 
the intermediate position the phase boundary separating 
it from PM is fully delineated below 70 kG. 

Given the clear in-plane anisotropy of the magneti- 
zation (fig. 5-10), it becomes desirable to systematically 
determine the angular dependence of the critical fields 
Hci.j(O) and locally saturated magnetizations Mj{9) of 
the HJlc metamagnetic transitions at T — 2 K. 

Fig. 11 shows a series of magnetization isotherms 

(T = 2 K) measured at various angles relative to the 
easy axis. The critical fields Hdj, for the transition be- 
tween states i and j, were determined from maxima in 
dM/dH, as exemplified in fig. 12, and are shown as full 
symbols in Fig. 13a. In most cases, Lorentzian fits of the 
derivative peaks were used (solid line in Fig. 12) to more 
accurately determine the critical field values. The open 
symbols in fig. 13 represent reflections of the measured 
data across the 9 = 0'^ direction; whenever the mea- 
sured points extend beyond the 0°..30° region, they al- 
most coincide with the calculated reflections, as expected 
for a symmetry direction. 

For 9 < 12° the antiferromagnetic AF ground state 
exists for fields up to about 20 kG, after which two 
closely spaced metamagnetic transitions occur, with crit- 
ical fields, at 61 = 0°, HcAfa = 20.5 kG (for the AF 
to Ml transition) and Hci,2 = 27.7 kG (corresponding 
to the transition from the Mi to the M2 state). A third 
transition from M2 to CL — SPM, around a critical field 
Hc2,CL-SPM = 53.7 kG, changes very little with angle 
up to ^ « 8°; for higher angles, another metamagnetic 
state M3 forms, being delineated by two distinct critical 
fields, Hc2.3 and Hc?,^cl-spm- As the former decreases 
with the angle, the latter soon reaches values around the 
maximum field of 55 kG available in the SQUID magne- 
tometer used for these measurements. In order to follow 
this latest transition in higher magnetic fields, additional 
measurements were taken in a different magnetometer, 
for fields up to 70 kG, and a slightly diflFerent tempera- 
ture (T = 1.85 K); these data are shown in the inset in 
Fig. 11a, but by 6* = 12°, Hc3,cl~spm becomes larger 
than 70 kG, therefore we can only anticipate that this 
transition still exists for larger angles. (The H — T phase 
diagrams (Fig. 9b and c) seem to indicate that this crit- 
ical field value increases from ~ 68 kG for 9 = 12°, to 
more than 140 kG at ^ = 30°). Slight differences can 
be noticed between the data sets taken in the two ma- 
chines, very likely due to the different temperatures at 



5 



which they were taken. A Hnear scahng of the two data 
sets by a factor of ^ 1.07 was necessary, for both the 
magnetization values and the critical fields; the scaling 
of the magnetization values can be explained by an as- 
sumption of slightly different angles between the applied 
field and the rotator axis in each magnetometer, while 
the field values may have changed with T according to 
the phase diagrams in Fig. 9. 

After the scaling of the two data sets, and after addi- 
tional calibration to the measurements on the large mass 
13.55 mg piece, the locally saturated magnetization val- 
ues were determined. The criterium used for determin- 
ing the magnetization for each state {Mj), was the onset 
M{H) value (fig. 12), i.e. the intersection of the linear 
fit of the Mj2 magnetization plateau and the highest- 
slope linear fit of the M{H) curve during the Mj-^ to 
Mj^ transition. More attention was given to determin- 
ing the magnetization for the first state (Mi), due to the 
limited field range over which this state exists. Several 
criteria tried in this case (onset value, midpoint between 
transitions, minimum in dM / dH or midpoint on the ap- 
propriate linear region on the M{H) curves) resulted in 
almost identical angular dependencies of Mi; moreover, 
using any of the aforementioned criteria, we were still 
unable to follow this state in the M{H) curves for angles 
beyond 25°. 

For 6 > 12°, the similar two sets of measurements are 
shown in Fig. lib. The same criteria were used for the de- 
termination of Mj (6) and H^j {9) . The two lower meta- 
magnetic transitions can also be seen in this region, while 
of the higher two, only -ffc2,3 is within our field range; as 
the angle increases, the first two transitions move closer 
in field {HcAF,i increases, while Hci,2 doesn't vary sig- 
nificantly with 9), such that the Mi state becomes very 
narrow, making its determination very difficult. M2 and 
M3 however appear as well defined plateaus, continuously 
decreasing, and increasing respectively, from the local ex- 
tremum values seen a.t 9 = 0°. Hc2.3 has a minimum of 
38.9 kG around 9 = 30°. 

The magnetization curves revealed four metamagnetic 
states, and their angular dependence is presented in 
Fig. 13b: Mi, M2 and Mcl-spm, which have local max- 
ima at 61 = 0° around 1.08 /lij/Tb, 3.00 ^s/Tb, and 
6.25 /is/Tb respectively, and M3 which exists only be- 
yond 9 — 8^ and has a maximum of 5.06 /is/Tb at 
9 = 30°. Similar to the Hdj in Fig. 13a, the open sym- 
bols in Fig. 13b represent reflections of the measured data 
across the 6* = 0° (easy axis) direction. 

The dotted lines in Fig. 13a and b are fits to Hcij{9) 
and Mj (9) respectively, as calculated based on the model 
that will be discussed below. Their angular dependencies 
are described by 1 / cos{9 ± ipj), and cos{9 ± ipj) re- 
spectively, with ipj — 0°, 30° or 60°. These values 
are integer or half-integer multiples of 360° / n, where 
n = 6 in our hexagonal system. Considering the six- 
fold symmetry of the this compound, these simple geo- 
metrical relationships render TbPtIn as very similar to 
RNi2B2C'^'"' or DyAgSbz'^, tetragonal compounds where 



the analogues (pj values were 0°, 45° or 90° (integer or 
half-integer multiples of 360° / n, where n = 4). 

As seen earlier in the cases of the tetragonal com- 
pounds HoNi2B2C^ or DyAgSbz^, simple angular depen- 
dencies of the critical fields, as well as of the locally sat- 
urated magnetizations exist in the hexagonal compound 
TbPtIn; this will be further confirmed by similar geomet- 
rical relationships that appear to exist in TmAgGe. 



IV. DATA ANALYSIS 

The field and temperature dependent magnetization 
measurements on TbPtIn (fig. 2 and 5) have shown that 
this compound is extremely anisotropic, with the mag- 
netic moments confined to the hexagonal basal plane. 
Moreover, when the direction of the applied field is var- 
ied within the basal plane, six fold anisotropy of the sat- 
urated magnetization is revealed in both TbPtIn and its 
dilution (Tb2;Yi_2:)PtIn (Fig. 6). Consequently, detailed 
magnetization measurements with i?_Lc were performed, 
allowing us to quantitatively describe the angular depen- 
dencies of the critical fields Hdj and locally saturated 
magnetizations Mj (Fig. 13). By analogy to the four- 
position clock model^i^ for tetragonal systems, we are 
now proposing a simple model for the net distribution of 
moments in the hexagonal compound TbPtIn: three co- 
planar Ising-like systems, 60° apart in the basal plane. 
Such a hypothesis was first suggested by the high field 
magnetization values observed in the pure compound 
TbPtIn, as weU as in the highly diluted (Tba;Yi_a;)PtIn. 
As the maximum measured magnetization for TbPtIn 
was around 6 /is /Tb^^^ (far smaller than the calculated 
9 /is value), it is reasonable to assume the existence of 
more metamagnetic transitions beyond our maximum ap- 
plied field H ~ 140 kC However, in the highly di- 
luted compound, where, within our field and tempera- 
ture ranges, we are only probing the paramagnetic state, 
the magnetization also reaches only ~ 6 /iB/Tb'^+ at the 
highest H. This is consistent with the M = 6 /iB/Tb'^+ 
corresponding to a crystal-field limited saturated para- 
magnetic (CL-SPM) state. Consequently we chose our 
model based on three Ising-systems such that it described 
the hexagonal symmetry of the compound having three 
magnetic ions in orthorhombic point symmetry, with the 
above value corresponding to saturation (in the limit of 
high CEF energy). In order to verify this hypothesis, 
the expected angular dependencies of possible moment 
configurations resulting from such a model will be com- 
pared with our measurements. Furthermore, we will use 
our experimental results to refine the model, by consid- 
ering multiples of the three Ising-like systems, resulting 
in more complex angular dependencies of the calculated 
magnetization and critical field values. 

In the P62m space group, TbPtIn assumes a hexag- 
onal crystal structure, with 3 Tb^^ ions at equivalent 
3g (orthorhombic) sites. The fact that a strong CEF 
anisotropy confines the local moments to the basal plane 
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calls for a two dimensional model, greatly simplifying 
the analysis. (A schematic description of an equiva- 
lent three-dimensional model has been introduced for 
DyAgGc''). Having three equivalent magnetic moments 
in orthorhombic point symmetry, one possible way to 
achieve the overall hexagonal symmetry is by restrict- 
ing the moments to three of the six-fold symmetry axes, 
60° apart, while allowing for both the 'up' (solid arrows) 
and 'down' (dotted arrows) positions for a given direc- 
tion (Fig. 14). Any specific Tb-site would, at low temper- 
atures, behave like an Ising system, with each third of the 
sites having parallel Ising directions. Thus, each mctam- 
agnetic state of TbPtIn can be described by a multiple 
S of three Ising-like systems along three [120] equivalent 
directions (the easy axes for this system). We will use \ 
, "f and / symbols, to denote the orientation of the three 
moments in their 'up' positions, and \ , ], and y sym- 
bols respectively, for the corresponding 'down' positions. 
The order of the arrows is not meaningful for our model; 
only the number of arrows for each orientation is signif- 
icant for the net distribution of moments. Moreover, we 
describe each metamagnetic state with the minimum-S 
value moment configuration consistent with the experi- 
mental data. However, higher S values are possible for 
most of the states, and information about the wave vec- 
tors (e.(7. from scattering experiments) would be required 
to determine unique S values. 

Since in our experiments we only measure the projec- 
tion of the magnetic moment along the field direction, the 
angular dependence of the magnetization Mj per moment 
of an arbitrary 3 S-moments configuration is 

Mj(0)///,,t(Tb3+) = [ Ef=i * cos(0 - 60°) + 

Sf=i *cosQ -f Yfi=\ ™'i * cos(6' -I- 60°) ] 

where 6* is a continuous variable representing the an- 
gle between the applied field and the closest easy axis 
(—30° < 9 < 30°), and the three sums give the mag- 
netization value due to each of the three directions of the 
Ising-like systems; the parameters equal ± 1, depend- 
ing on whether a certain moment is in the 'up' (+ 1) or 
'down' (— 1) position for the respective direction. We 
restrict our model description to the 0° < < 30° 
angular region, which, by symmetry across the — 0° 
direction, also describes the —30° < 6* < 0° region. 

For 5 = 1, our model corresponds to one set of three 
such Ising-like systems. We assume that in high applied 
fields, the three magnetic moments occupy the three al- 
lowed easy axes closest to the direction of the field; as 
the field is lowered, the metamagnetic transitions oc- 
cur such that the measured magnetization is being de- 
creased with H. In this hypothesis, there are three dis- 
tinct moment configurations for the system: (\t /*) for 
the CL-SPM state, (\t/') for intermediate field values, 
and (\Ti/) for the AF ground state. The above formula 
yields the following angular dependencies of the result- 
ing longitudinal (measured) magnetizations: 2/3* cos 6*, 
2/3*cos(0 — 60°) and respectively, represented by open 
circles in Fig. 15. 

Fig. 13b shows that such a model only describes the 



CL-SPM state McL-sPAf(6')//^sat(Tb3+) = 2/3*cos6' 
and the AF ground state Maf / pi-sati^h^^) = of 
TbPtIn; according to the proposed model, the local mo- 
ment configurations from which the above angular de- 
pendencies follow are (\t y) and (\ti/) respectively. 
It is worth noting that the CL-SPM magnetization value, 
calculated based on the above moment configuration, is 
6 fj,B, smaller than the measured 6.25 /is- One possi- 
ble explanation is that with increasing field, the system 
is slowly approaching the CEF splitting energy. This is 
also consistent with the increasing plateaus in the high 
field magnetization data in Fig. 5a; however, the extrapo- 
lation of these plateaus down to iJ = results in smaller 
values for the [120] direction (6.13 (Ib) and the [110] di- 
rection (5.35 /is), closer to the calculated values (6.0 /is 
and \/3/2 * 6 /is — 5.2 fis respectively). 

To characterize all the other observed metamagnetic 
states, larger S-values are needed, i.e. the local mo- 
ment configurations are described by an integer multiple 
S > 1 of sets of three Ising-like systems. Fig. 15 (crosses) 
also shows all possible angular dependencies of the mag- 
netizations resulting from such a generalized model when 
S — 2. By comparison with the experimental data, it 
appears from Fig. 13b that two more metamagnetic states 
can now be described with S = 2: 

M3(0)//isat(Tb3+) = 2V3/6 * cos(6l - 30°) (moment 
configuration (\\TT//) ), 

and 

M2{0)/^sat{Th'^^) — 2/6 * COS (moment configura- 
tion (WTT//) ) . 

There is still one more metamagnetic state, AIi, 
which cannot be described within the 5 = 2 model; 
however, for most of the angular range, its magne- 
tization has an angular dependence consistent with: 
Mi(6')//i5at(Tb3+) = 2/18*cos6'. The cose* depen- 
dence (i.e. cos is an even function) requires that the 
moment configuration be symmetric with respect to the 

= 0° direction. The simplest possibility is (\t/'), 
for which the magnetization varies as 2/3 * cos 9; to this, 
a number of sets of three moments needs to be added, 
with zero net magnetization (e.g. multiples of (\Ti/) 
or (Wtiv//')), to get a resulting magnetization 
amplitude of 2/18. Consequently a minimum 5 = 6 
configuration (\\\\\\TTTTTT//////) or 
(WWWiiTTTT//////) yields the desired 
calculated magnetization 2/18 * cos 9. 

Assuming the above net distributions of moments for 
the observed metamagnetic states, one can derive the ex- 
pected angular dependencies of the critical fields. Com- 
parison between the data in Fig. 13a and these calculated 
Hci,j (9) values will further confirm or refute the net dis- 
tributions of moments proposed above. 

Since the energy associated with a magnetic moment 
M in an applied field H is M ■ , the corresponding 
energy difference AEji between metamagnetic states Mi 
and Mj is: 

AEji = NIj - ll -M,-H 
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If there is a critical energy Ec = AEji to be exceeded 
for a metamagnetic transition between states i and j to 
occur, than the critical field value is given by: 

rr _ Eg 

where Mj and Mi are the measured (projections along 
the field) respective magnetizations. The numerator in 
the above expression is angle and field independent, and 
the angular dependence of Hdj follows only from the 
denominator. In other words, 

rr 1 

^cij- Mj-M,- 

Consequently, the expected critical field values, shown 
as dotted hues in Fig. 13a, are: 

HcAF.m - 1/COS0, 

HciM^) ^ 1/cos^, 
Hc2,m ^ l/cos(0-6OO), 
Hc3,CL-SPM {6) - 1/008(6* + 60") 
and Hc2,cl-spm{6) ^ 1/ co&O. 

The reflections across the 9 — 0^ direction result 
from the above formulas, when substituting 9 with —9] 
moreover, since cos is an even function, this is equivalent 
to a change in sign only for the Lp in the above expressions 
written as 1/ cos{9 — ip). 

As described above and similar to the analogous study 
in the tetragonal compound HoNi2B2C^, in most cases 
we used maxima in dM/dH, and not the on-set criterion, 
to determine the critical field values, because the mag- 
netizations during the transition were not always linear; 
however, comparison with the calculated critical fields 
based on the above model is still appropriate, given that 
only small departures from linearity were encountered, 
mostly close to the bordering states {Mi and Mj). (The 
non-linear change of the magnetization with increasing 
field indicates that other factors {i.e. the demagnetiza- 
tion factor of the sampleii, coexistence of more than two 
phases, non-linear superposition of the various states) 
may be responsible for the broadening of the transition) . 

Comparison of the measured critical fields and lo- 
cally saturated magnetizations (Fig. 13a, b full and open 
symbols) with the calculated values as described above 
(Fig. 13a, b dotted lines) confirms, in most cases, the 
assumed local distribution of moments. However, the 
first metamagnetic state Mi follows well the calculated 
2/18*Aisat(Tb3+)*cos(6') dependence up to 6* « 25°, after 
which it is difficult to determine it with reasonable accu- 
racy. Despite the fact that both HcAf.i and i?ci,2 should 
depend on this magnetization value, the former follows 
the expected angular dependence fairly well, whereas the 
latter falls under the calculated l/cos0 curve. Even 
though at this point we don't have a rigorous calculation 
to support our assumption, we anticipate that the prox- 
imity of the two lowest transitions requires a calculation 
with more than two coexisting phases, which may render 
a better fit of the observed experimental data in this re- 
gion. On the other hand, the i?c2,3 and Hcs^cl-spm crit- 
ical fields, which have the most evident angular depen- 
dence, are well fitted by the calculated functions based 
on the present model. It is thus reasonable to assume 
that, whereas possibilities for refining the model exist. 



in the simple form that we present here it describes our 
system fairly well. 

The polar plot in Fig. 16 helps in understanding how 
the hexagonal crystal structure of this compound is re- 
fiected in the angular dependence of the metamagnetic 
phase transitions: similar to the polar phase diagrams for 
tetragonal compounds IIoNi2B2C— or DyAgSb2^, when 
we plot He * sin 9 vs. He * cos 9, the phase boundaries be- 
come straight lines, with slopes equal to either ± l/v^ or 
oo. These slopes correspond to directions either parallel 
or perpendicular to the high symmetry axes {i.e. [110] or 
[120]) in the hexagonal structure, just as within the four 
position clock model the corresponding phase lines were 
either parallel or perpendicular to the tetragonal high 
symmetry axes ([110] or [010]). As described in Myers 
et al."^ , the equations of these straight lines in polar co- 
ordinates can be used to verify the transitions already 
discussed: if we substitute the above slope values in the 
general formula 

R{9) = a/(sin6'-6*cos6i), 

for a line with slope b, we get: 

Hc{9) - l/cos6l for 6 = oo 

Ha{9) - l/cos(6'- 60") for = - l/\/3 

or 

He{9) - 1/008(6* + 60") for 6 = 1/^/3, 
which are consistent with the angular dependencies of 
the transitions determined above. 

For the most part, the experimental points fall onto 
the calculated straight lines, as expected. Some devi- 
ations from the straight lines can be noticed, with the 
most evident one for Hci.2, for which we already empha- 
sized the necessity of a more complex model. In a similar 
manner (even though only for fewer angles), -ffc2.3 curves 
under the calculated straight line as we move away from 
the [120] easy axis, while some even smaller deviations 
from linearity are apparent in HcAF,i] ^^^^ may indicate 
that special attention needs to be paid in determining the 
angular dependence of the critical fields, when rotating 
from the proximity of one easy axis to another. 



V. TmAgGe 

We already reported the basic magnetic properties of 
TmAgG&£, which strongly resemble those of TbPtIn: the 
magnetic susceptibility is extremely anisotropic (Fig. 17), 
indicating antiferromagnetic order below Tn = 4.2 K. 
The local magnetic moments are confined by the strong 
CEF anisotropy to the basal plane, both below and above 
Tn. This can also be seen in the field dependent magne- 
tization measurements, shown in Fig. 18, where in-plane 
anisotropy of the ordered state is also apparent. Similar 
to the case of TbPtIn, several metamagnetic transitions 
exist for both i?||[110] and i/||[120]; these result in mag- 
netization values of 4.92 and 4.30 /is respectively, at 
H ^ 70 kC, far below ^sat(Tni'^+) = 7.0 /is, whereas 
for the c direction, the magnetization is linear and much 
smaller up to the maximum applied field. The ratio of the 
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two in-plane magnetizations is J^^^^^-^Qj^ = 0.87, close to 

the cos 30° value expected within the model described be- 
fore for TbPtIn for the CL-SPM state. However, the two 
absolute values are larger than the corresponding ones, 
calculated from the above model: M([110]) = A.&T^b 
and M([120]) = 4.00 hb, but the extrapolation of the 
high- field plateaus down to = (solid lines in Fig. 18) 
yields magnetizations very close to these calculated val- 
ues. As in the case of TbPtIn, the slight increase of 
the magnetization plateaus after the supposed saturation 
may be caused by the slow approach of the CEF splitting 
energy. 

In calculating the above expected magnetization val- 
ues, we assumed the easy axes to be along the [110]- 
cquivalcnt directions, based on the directions where max- 
imum magnetization values at if = 70 kG were achieved 
(fig. 18). This is consistent with the angular dependent 
magnetization mcasiircmcnt for H = 70 kG shown in 
Fig. 19, where angle 9 was measured from the [110] di- 
rection; thus the six-fold symmetric magnetization has 
maxima occurring for the [110]-cquivalcnt directions {i.e. 
for 6 = n * GO'^, where n is an integer). Consequently, 
for TmAgGe the angle 9 will be measured from the clos- 
est [110] easy axis. The comparison of TmAgGe and 
TbPtIn indicates that, even though the easy axes in the 
two compounds correspond to the two different sets of 
six- fold symmetry directions, as we shall see, their phys- 
ical properties are very similar. 

Around each easy axis, these magnetization measure- 
ments follow the cos 9 angular dependence (the solid 
line in Fig. 19), as expected within our proposed model. 
Some differences between the experimental data (filled 
circles in fig. 19) and the calculated magnetization could 
be caused by small misalignment of the sample (ren- 
dering slightly asymmetric measured peaks), or by the 
strong interactions between the local moments. Similar 
to TbPtIn, this also indicates that different metamag- 
netic states cross H = 70 kG at different angles. 

From the M(T) \h,0 (fig.20) and M{H) \T,e (fig-21) 
mciasurcments, detailed H ~ T phase diagrams for this 
compound can be determined. They are shown in 
Fig.22a,c, for field along the [110] or 6' = 0°, and [120] or 
9 ~ 30*^ directions respectively, with an intermediate- 
position 9 ~ 24P phase diagram in Fig. 22b. In the 
same manner used for TbPtIn, the points in these phase 
diagrams have been obtained from maxima in either 
d{M * T/H)/dT for constant field (fuU symbols) or in 
dM/dH for fixed temperatures (open symbols). 

For if||[110] (fig. 22a), at low temperatures the antifer- 
romagnetic, AF, ground state exists for H < 3.1 kG, 
followed by a small intermediate phase Mi (up to 4.4 
kG) and a larger state M2 above. This latest phase ex- 
tends up to 8.9 kG, after which, at low temperatures, the 
system reaches the crystal field-limited saturated param- 
agnetic, CL-SPM, state. As temperature is increased, 
the Ml phase disappears around 2.5 K , and a direct 
transition from the AF to the M2 state occurs at a de- 
creasing critical field value. The upper phase boundary 



(for the M2 to the PM state transition) also falls down 
in field as T increases, such that at very low fields only 
one transition is observed close to Tj^ = 4.2 K. 

As we rotate away from the easy axis, the low-field 
phase diagram changes very little, with a small enhance- 
ment of the critical field values towards low temperatures. 
As field is being increased, M2 is getting smaller as a 
new distinct phase M3 forms. Its upper bordering line 
appears to have a strong angular dependence, as can be 
seen in Fig. 22b and c, similar to the upper-most phase 
boundary seen in TbPtIn : for i?||[120] (Fig.22c), this 
phase boundary is an almost vertical line at T « 4.5K, 
up to our maximum applied field H = 70 kG. As a 
consequence, the corresponding points on this line have 
been determined from d{M * T/H)/dT data, as shown 
in Fig. 20c for high fields, and could not be identified in 
the field-dependent derivatives. At low temperatures, the 
HcAF,i and ifci,2 values (3.62 kG and 4.86 kG respec- 
tively) are very close to the corresponding ones in the 
[110] direction, whereas the M2 to M3 transition occurs 
around 7.0 kG. These three phase lines merge around 
T = 3.0 K, such that for higher temperatures a single 
transition occurs at decreasing fields. This line appears 
to intersect the H = axis around Tjv = 4.2 K. 

The intermediate-orientation phase diagram presented 
in Fig. 22b allowed us to observe the upper- most phase 
line moving down in field at low temperatures, such that 
for 9 Ri 24°, it intersects the T = axis close to 
40.1 kG. In this orientation, this phase boimdary can be 
identified in the M(H) derivative, as shown in Fig. 23 for 
T = 1.85 K. However, the high field peak in dM/dH 
is poorly defined, making the determination of the corre- 
sponding critical field value more difficult. For a more 
precise estimate, another criterium was used together 
with the derivative maxima, as illustrated in the inset 
in Fig. 23: the mid-point (large dot) on the highest-slope 
linear fit (solid line) of the magnetization data around 
the transition. Also shown is the error bar for this criti- 
cal field value, as determined from the two criteria used 
here. 

The fact that this line is now apparent in both M{T) 
and M{H) data is further confirmation that this phase 
boundary exists, whereas at lower fields, the only no- 
ticeable difference from the i?||[120] direction is the per- 
sistence of the M2 state up to higher {i.e. ~ 3.5 K) 
temperatures. 

A number of similarities between TbPtIn and TmAgGe 
have already been established: same crystal structure, 
antiferromagnetic ground state, extremely anisotropic 
magnetization, in-plane anisotropy and metamagnetism 
leading to crystal field-limited saturated magnetizations 
smaller than the calculated single ion Usat values. As 
a consequence, we proceed to study the angular depen- 
dence of the planar metamagnetism in TmAgGe, and 
subsequently apply the model developed for TbPtIn to 
the case of this compound. 

When we fix the temperature at T = 2 K, the angu- 
lar dependence of the metamagnetic transitions can be 
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studied based on the M{H) isotherms shown in Fig. 24. 
The critical fields and the locally saturated magnetiza- 
tion values (full symbols in fig. 25) have been determined 
as maxima in dM/dH, and from on-set values respec- 
tively (see the TbPtIn section). An exception was made 
for McL-SPM above 10°, and the criterion used for de- 
termining this state is described below. Moreover, be- 
cause of the proximity of the first two transitions, the 
Ml state is poorly defined; no precise saturated magne- 
tization data could be extracted for this state, but the 
phase diagrams in Fig. 22, as well as the angular depen- 
dent critical fields in Fig. 25a, arc consistent with the ex- 
istence of this phase. As before, the open symbols in 
Fig. 25 represent reflections of the measured data across 
the 6 = 0° direction. The TmAgGe measurements 
allowed us to determine the critical field and locally sat- 
urated magnetization values for the full angular range 
(—30° < < 30°); the resulting somewhat asymmet- 
ric data (most obvious in the case of the M2 data) maybe 
due to a small sample misalignment. The experimental 
data (full symbols) together with the reflections (open 
symbols) in Fig. 25 give the caliper of the error bars for 
these measurements. 

For 6 < 10°, two closely spaced metamagnetic 
transitions can be seen in Fig.25a, with critical fields 
HcAF.i = 3.0 kG and -ffci.2 = 4.37 kG respectively at 
^ = 0°, followed by a third transition Hc2,cl-spm at 
~ 9.37 kG. It should be noted that these i?ci ,j values are 
slightly different from the corresponding ones (3.10 kG, 
4.36 kG and 8.92 kG respectively at T = 2.0 K, 6* = 0°) 
in the H~T phase diagrams (Fig. 22a), as they have been 
determined from two distinct measurements. Thus small 
errors in the angular position (± 1°) may convert into 
small errors in the critical field values (< 3%). 

Somewhat larger differences between the two data sets 
are observed for 6 > 10°, specifically for Hc3,cl-spm, 
which varies more rapidly with the angle than any other 
critical field. In this angular region, the two lower tran- 
sitions occur at almost the same critical fields as below 
10°, whereas the critical field for the third one slowly 
decreases with angle, as a fourth transition appears and 
rapidly moves up in field. Consequently, the Mi metam- 
agnetic state changes little with the angle, whereas the 
M2 state narrows down as the bordering critical fields 
move closer to each other. The fourth transition being 
very broad makes the determination of the M3 state fairly 
difficult. Also, with Hc3.cl-spm broadening out and 
rapidly moving towards our field limit {i.e. 70 kG), it 
was difficult to get a meaningful linear fit of the M{H) 
curves during the M3 to CL — SPM transition; instead 
wc used the intersection of the maximum slope line cor- 
responding to the M2 to M3 transition, and the best lin- 
ear fit of the highest magnetization state, to determine 
McL-SPM for e > 10° (Fig.24b, inset). 

The best fits to the experimentally measured angular 
dependent data are shown in Fig. 25 as dotted lines. We 
will use these fits to infer the net distribution of moments 
as multiples S of three Ising-like systems, similar to the 



case of TbPtln. As mentioned before, we infer that the 
Ml state should exist based on the angular dependent 
critical fields in Fig. 25a, and the T = 2.0 K metam- 
agnetic phases revealed by the phase diagrams in Fig. 22. 
Consequently, in Fig. 25b we are only showing the ex- 
pected angular dependence of such a phase, by analogy 
with the TbPtIn case: Mf"'7/isat(Tm3+) = 2/18*cos0, 
which appears to be the upper limit of these magneti- 
zation values, as indicated by the error bars shown in 
Fig. 25b. As already seen for TbPtIn, the moment con- 
figuration that would result in such a Mi magnetization 
isaS = estate: (\\\\\\TTTTTT//////) or 
(WWWiiTTTT//////). A'hiO) has a max- 
imum value at = 0° equal to 2.37 /iB/Tm'^+, close 
to 2/6 * Hsat{Tm^+) = 2/6 * 7/iB/Tm3+, and a cos 61 
angular dependence. This suggests that a possible net 
distribution of moments for this state, realized with a 
minimum S = 2, could be (WTT,//')- If one local 
moment is flipped from its ^ position to Z' in the pre- 
vious state, the resulting state could be described by the 
(WTT^^) configuration, whose magnetization varies 
as 2\/3/6*/i,,at(Tm3+)*cos(6'-30°); this fits wefl the mea- 
sured M3 data, indicating that the previously assumed 
local moment distribution may be appropriate for this 
metamagnetic state. When all the magnetic moments 
are in their 'up' positions, the CL-SPM (WTT/'/') 
state is achieved, and the corresponding angular depen- 
dence is 4/6 * /Ltsat(Tm^+) * cos9. This best describes 
the last observed metamagnetic state, which has a max- 
imum of 4.68y^B/Tm^+ around 6 = 0°, very close to 
4/6* jUsat(Tm3+) = 4/6 * 7/ZB/Tm3+. (As before, the 
order of the arrows used to describe the net distribution 
of moments has no physical meaning). 

According to the calculation given in the case of 
TbPtIn, for TmAgGe one would also expect the critical 
fields to vary with the angle 9 as 



Using the net distributions of moments assumed above 
to best describe the locally saturated magnetization 
states, we expect the following angular dependencies of 
the critical fields: 

HcAfAO) ~ 1/cos^, 

Hci.2{0) ~ l/cos6' 

He2,m - 1/COS(0-6O°), 
Hc3,CL-SPM ^ l/cos(6l-h60°) 

and Hc2,cl-spm{9) ~ l/cosff. 

The experimental data (full symbols in Fig. 25a) are 
in good agreement with these calculated critical fields, 
with Hc2,cl-spm{6) present only for 6 < 10°, while 
the i?c2,3 and Hcz^cl-spm exist only for 6 > 10°. 
This is consistent with the presence of the M3 state for 
angles larger than 10°, even though experimentally we 
were only able to accurately determine it for 9 > 16°. 
(The bordering transitions of this state are very close in 
field when 9 is close to 10°, which made the determination 
of M3 in this angular region difficult). It is worth pointing 
out the excellent fit of the measured Hcz.cl-spm data 
with the calculated angular dependency, this transition 
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showing the most dramatic change with the angle 9. 

Apart from the absolute values of the critical fields 
and the locally saturated magnetizations, the Hc{9) and 
^sat{0) phase diagrams for TmAgGe (Fig. 25) are identi- 
cal to the TbPtIn analogues in Fig. 13: to the same num- 
ber of critical fields with identical angular dependencies 
correspond identical metamagnetic values (scaled to the 
saturated moment of the respective R^"*" ion), which also 
vary similarly with the angle. This is consistent with our 
model being indeed a general description of the Fe2P- 
type systems, or even more generally, of hexagonal sys- 
tems with the R in orthorhombic point symmetry. 

When the Hc{6) phase diagram for TmAgGe is con- 
verted into a polar plot (Fig. 26), similar to that for 
TbPtIn, we again notice that the phase boundaries are 
straight lines, with ± l/>/3 or oo slopes. From the equa- 
tions of these straight lines, we can once more confirm 
the transitions determined before: 

Hcie) ~ l/cos6lfor6 = oo 

H^le) ~ l/cosC^i-eoO) for 6 = - l/\/3 

or 

Hcie) ~ l/cos(6'-F600) for & = l/\/3. 

As already noted, in the case of TmAgGe, we were 
able to determine the critical fields from experimental 
data for ^ = — 30'^ .. 30°, as seen in the phase diagrams 
in Fig. 25a, as well as in Fig. 26. Slight differences between 
the expected straight lines in Fig. 26 and measured critical 
fields can be noticed for iJc2,3 for angles close to ± 30*^, 
or for Hc3,cL-SPM also for large angles; besides being a 
consequence of small misorientation of the sample, this 
may indicate, similar to the TbPtIn case, that a more 
complex model needs to be used to describe the regions 
around the 'hard' in-plane direction. 



VI. ANGULAR DEPENDENT 
METAMAGNETISM AT T = 20 K IN TbPtIn 

Whereas TmAgGe has fairly simple H — T phase di- 
agrams, with all metamagnetic phases present at low 
temperatures, the TbPtIn phase diagrams are somewhat 
more complex, manifesting an additional intermediate- 
temperature phase, M4 (Fig. 9). In order to perform a 
similar angular dependent study of this metamagnetic 
state, a set of M{H) \e data was taken at T = 20 K. 
Assuming that all existing transitions have been identi- 
fied and are shown in fig. 9 (at least below 20 K), at this 
temperature the M{H) curves should intersect the same 
magnetic phases as in the low temperature case, with the 
exception of Mi; instead, the measurements at T = 20 
K intersect the bubble-like phase M4, as seen in the three 
different orientations phase diagrams in Fig. 9. 

Fig.27 shows the M{H) isotherms at T = 20 K for 
various angles 9. The Mj(9) and Hci.j{9) phase diagrams 
have been determined as described before for TbPtIn or 
TmAgGe for T = 2 K, and are shown in fig.28. It 
should be noted that, due to the enhanced temperature, 
all transitions are broadened, and the locally saturated 



magnetization plateaus are no longer horizontal. Both 
of these facts make the analysis of these data somewhat 
harder and more ambiguous. 

For 9 < 12" (Fig. 28a), the lowest transition changes 
very little with angle, having a critical field value 
HcAF,4: ~ 8 kG. As field is being increased, two more 
transitions can be observed for angles lower than 8°, 
with local minima of the critical fields, at ^ = 0°, of 
Hc4,2 = 16.2 kG and Hc2,cl-spm = 48.7 kG re- 
spectively. For larger angles, the highest transition splits 
into two different ones, -ffc2,3, with decreasing values as 
we rotate away from the easy axis, and Hc3,cl-spm, 
which rapidly increases above our field limit (i.e. 55 kG) 
around 6 = 12^. It should be noted that we are still 
referring to the high-field state at T = 20 K as the crys- 
tal field-limited saturated paramagnetic CL-SPM state, 
even though it is possible that cross-over to the param- 
agnetic PM state has occurred between 2 K and 20 K at 
high H. (This would be a plausible explanation for the 
measured magnetization values for this high field state 
being, as seen below, lower than the calculated values.) 

The locally saturated magnetization of the M4 state 
is equal to ~ 0.25 /zs/Tb at 6* — 0° and doesn't ap- 
pear to change much with the angle. According to the 
H — T phase diagrams in Fig. 9, all higher metamagnetic 
states are identical at low (T = 2 K) temperature and 
at T = 20 K; consequently, they seem to have similar 
angular dependencies (fig.28b): M2 and Mcl-spm have 
local maxima around 3.3/iB/Tb and 5.7/iB/Tb respec- 
tively, at 6* = 0", and slowly decrease with increasing 
angle in this region. Beyond ^ 8" , a third metamagnetic 
state should exist, defined by the ifc2,3 and Hc3,cl-spm 
critical fields; however, it is difficult to identify it in this 
angular region, given the broadness of the bordering tran- 
sitions and the Hc3,cl-spm proximity to our field limit. 

As we move further away from the easy axis {i.e. 
6 > 12°, Fig. 28b), the lower two metamagnetic states 
can again be observed, whereas the CL-SPM state may 
still exist for fields larger than 55 kG (also apparent from 
the H — T phase diagram in Fig. 9b). Also, we can now 
see the third state M3 slowly increasing with angle 9, 
similar to the low temperature case. 

The T = 20 K phase diagrams (fig.28) are very sim- 
ilar to their low temperature analogues (Fig. 13), except 
for the M4 state, and some evident differences between 
the experimental data and the model calculations (dot- 
ted lines). It was rather difficult to determine the angular 
dependencies of M4 and M2, therefore we had to infer the 
possible theoretical fits in a more indirect way: as we al- 
ready mentioned, the state described by M2 at T = 20 K 
should be the same as the corresponding one at low tem- 
perature, since they characterize the same metamagnetic 
phase. Therefore we expect it to vary with the angle like 
2/6 * Msat(Tb3+) ^ pogg/ (fig.i3b). For large angles, this 
is consistent with the measured data in Fig. 28b, whereas 
significant deviations can be noticed closer to ^ = 0°. 
The lower magnetization state M4 has much smaller val- 
ues than any of the states characterized at T = 2 K, 
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therefore we cannot fit it accurately with a calculated an- 
gular dependence; however, i?c4,2 can be fitted with 16.2 
kG *1/ cos^ and should relate to M4 through 

HciA^) ~ i/[M2{e)-M4{e)] 

or 

l/cos6l - l/[2/6*/i^at(Tb3+)*cos6'-Af4(6')]. 

Thus M4{9) should vary like Msat,4 * cos 6*, with a lo- 
cally saturated magnetization MsatA ~ 0.25* fisat 
but, as already mentioned, it is difhcult to determine it 
with reasonable accuracy. However, the corresponding 
local moment configuration should be similar to the low 
temperature one, in order to get the cos9 dependence, ex- 
cept that the number S of three Ising-like systems which 
would yield the appropriate Msa,t,4 value is uncertain. 

All other locally saturated magnetizations and critical 
fields can be best fitted with the same angular dependen- 
cies as for the low temperature case: 

MsiO) = 2V3/6*//,at(Tb3+)*cos(0-3O°), 

McL-SPAliO) = 4/6*//,at(Tb3+)*COS0, 

and 

HcAFi ^ 1/ cos9, 
H,2,3 - l/cos(0-6O"), 

Hc2.CL-SPM ~ 1/ cos 6*, 

l/cos(6i-t-600). 

Apart from the already mentioned differences between 
the measured data and the calculated curves, small de- 
partures from the corresponding theoretical angular de- 
pendencies can be noticed for M^IO) and Hc2,3] a more 
significant difference appears for the saturated magneti- 
zation state Mcl-spm(Q), which seems to have the ex- 
pected angular dependence, but with smaller values than 
the calculated ones. This may be a high-temperature ef- 
fect {i.e. cross-over from low-T CL-SPM state to high-T 
paramagnetic PM state), or it may be one more indica- 
tion that a more refined model is needed. 

A polar plot analogues to the low temperature case 
(Fig.29) shows that at T = 20 K, the critical fields 
are still well described by straight lines, but with more 
pronounced differences between experiment and the the- 
oretical calculations. We attribute these to the thermal 
broadening at this temperature, but, as already seen in 
the low temperature requirement for a more com- 

plex model cannot be excluded. 



VII. SUMMARY 

Motivated by the extensive work done on highly 
anisotropic local moment systems with tetragonal unit 
cells and unique rare earth sites of tetragonal point 
symmetry, we have performed detailed studies on two 
highly anisotropic, local moment, hexagonal compounds: 
TbPtIn and TmAgGe. Whereas both of these com- 
pounds are ternary members of the Fe2P class of materi- 
als, they have different ligands. In addition, whereas both 
of these compounds manifest extreme planar anisotropy, 
they have different easy axes: [120] for TbPtIn and [110] 
for TmAgGe. Even with these differences we have found 



that these two compounds have very similar H — T as 
well as H — phase diagrams. 

TbPtIn and TmAgGe have a single rare earth site with 
orthorhombic point symmetry (with three rare earth sites 
per unit cell) and both compounds have high field satu- 
rated moments well below the single ion values. These 
two observations, combined with our experience with 
the four position clock model that was developed for 
the tetragonal compounds with rare earths in tetragonal 
symmetry, lead us to propose a similar model for these 
Fe2P-type compounds: a triple coplanar Ising model, 
which consists of three Ising-like moments per unit cell, 
with their Ising axes within the basal plane and rotated 
by 60" with respect to each other. This model preserves 
the six-fold symmetry at high fields and also explains 
why the saturated moments are significantly lower than 
the free ion values. By analyzing the magnitudes and an- 
gular dependencies of the critical metamagnetic fields, as 
well as the locally saturated magnetizations within the 
framework of this model, we can infer the net distribu- 
tion of moments along the six possible moment orien- 
tations. However, field-dependent neutron diffraction or 
magnetic x-ray measurements are needed to test these 
hypothetical net distributions of moments and to obtain 
the propagation vectors for each magnetically ordered 
state. (Some preliminary experiments on TbPtIn are 
already reportediS:, confirming the two low-field order- 
ing temperatures that we observed in the magnetization 
measurements. Moreover, the neutron data are consis- 
tent with our inferred direction of the moments.) 

The successful extension of the four state clock model 
to the current triple coplanar Ising-like model implies 
that a wider set of local moment compounds with planar 
anisotropy can be understood in a similarly simple man- 
ner. Clearly tetragonal unit cell compounds with the rare 
earth in orthorhombic point symmetry could be expected 
to behave in a manner similar to TbPtIn and TmAgGe, 
i.e. to form a class of double coplanar Ising model mate- 
rials. In a similar manner, hexagonal unit cell compounds 
with the rare earth in hexagonal point symmetry could 
be expected to behave in a manner similar to HoNi2B2C^ 
or DyAgSb^^, i.e. to form a class of six position clock 
model materials. Hexagonal unit cell materials may of- 
fer one other, potentially new class of materials: highly 
planar anisotropic compounds with the rare earth in trig- 
onal point symmetry. In this case we anticipate a double 
coplanar three position clock model. Such compounds 
would have two magnetic sites, each with three possible 
positions 120" apart. Fig. 30 gives schematic representa- 
tions of all expected five models described above. 
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FIG. 1: Projection of the ZrNiAl-type crystal structure along 
the hexagonal c-axis: R -large circles, M (Pt or Ag)-medium 
circles, X (In or Ge)- small circles. Light circles: z = 
plane, dark circles: z = 1/2 plane. (Note: the Pt or Ag 
atoms (medium circles) appear in both planes.) 




FIG. 2: Anisotropic inverse susceptibility of TbPtIn (sym- 
bols) and the calculated average (line); inset: low- 
temperature anisotropic susceptibilities. 
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FIG. 3: (a) Low-temperature d{xave *T)/dT, with vertical 
dotted lines marking the peaks positions; (b) specific heat 
Cp(T); inset: magnetic entropy SmiT); (c) low-temperature 
resistivity p(T), for current flowing in the basal plane {i \\ ab). 
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FIG. 4: Anisotropic inverse susceptibility (symbols) and cal- 
culated average (line) for H = 1 kG; upper inset: low- temper- 
ature anisotropic susceptibility; lower inset: field-dependent 
anisotropic magnetization, for T = 2 K. 
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FIG. 5: (a) Anisotropic field-dependent magnetization (in- 
creasing and decreasing field data indicated by arrows) and 
(b) transverse magnetoresistance, for increasing (symbols) 
and decreasing (line) field. 
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FIG. 6: M(6') of TbPtIn (full circles) and Tbo.o25Yo.975PtIn 
(open circles) for T = 2 K and H = 55 kG, H_Lc. Solid line 
represents the calculated Mmax * cos{0 — n* 60°), n-integer. 
Note: = O" is defined at the [120] direction. 
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FIG. 7; M(T) data for various fields for two in-plane orien- 
tations of the applied field: (a) H||[120] for H = 1, 2, 7.5, 15, 
20, 25, 30, 40, 45, 50 and 55 kG, and (b) H||[110] for H = 
1, 2.5, 3.5, 5, 7.5, 10, 15, 20, 26, 30, 37.5, 40, 42.5, 45, 50 
and 55 kG. Insets show enlarged M(T)*T derivatives (dotted 
lines) for H = 20 kG, together with the Lorentzian fits of the 
maxima (solid lines) , to cxamplify how the points represented 
by full symbols on the H-T phase diagrams were determined. 
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FIG. 8: M(H) isotherms for two in-plaric orientations of the 
apphed field: (a) H||[120] for T = 3, 5, 10, 15, 17.5, 20, 25, 
30, 35, 37.5, 40 and 45 K, and (b) H||[110] for T = 2, 5, 
10, 15, 20, 25, 30, 35, 40, 42.5, 45 and 47.5 K. Insets show 
enlarged M(H) derivatives (dotted lines) for T = 20 K and 
30 K, for the corresponding field directions, together with the 
Lorentzian fits of the maxima (solid linos), to oxamplify how 
the points represented by open symbols on the H-T phase 
diagrams were determined. 



21 



70 



60 



50 



40 



30 



20 



10 



(a) 



• from M(T) data 
O from M(H) data 



H||[120] 



PM 



M. 



AF 



o 



oo 



O ^ AF" 




M, 
AF 



On O 
o°»ocpoO 



AF' • 



(c) 



H||[110] 




^3 

8° o • o 

M. 



J o O » o 
M 

i* o 



• \ 

"-in* 



10 20 30 40 50 



M, 



O Oo 

Oot) 



PM 



AF 



o 



o 



O 



10 



20 30 
T(K) 



40 



50 



10 



20 30 
T(K) 



40 



50 



10 



20 



30 
T(K) 



40 



50 



60 



FIG. 9: H-T phase diagrams for (a) H||[120] and (c) H||[110], 
as determined from the magnetization data in Fig.7-8, with an 
intermediate field orientation phase diagram (see text) shown 
in (b). 




FIG. 10: (a) M(T) data for various fields: H = 1 kG and 
5-70 kG (AH = 5 kG) for the intermediate field direction (see 
text). The small arrows indicate the highest-field transitions 
at the corresponding temperature, (b) M(H) for T = 10 K, 
and the corresponding derivative, as an example of how the 
open symbols on the upper-most phase line in Fig. 9b have 
been determined. 
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FIG. 11: M(H) isotherms at T = 2.0 K for (a) < 61 < 12" 
and (b) 12° < 9 < 30° (A6' = 1°); inset: enlarged high 
field, T = 1.85 K (see text) isotherms. Arrows indicate the 
direction of increasing 9. 
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FIG. 12: Enlarged M(H; 9 = 12") plot, and the correspond- 
ing derivative, illustrating the criteria used in determining 
the points in Fig. 13: large dot on the dM/dH plot indicates 
the critical field value, determined by the maximum of the 
Lorentzian fit (solid line) of the peak; straight lines on the 
M(H) plot are fits to the magnetization plateau (thin), ex- 
trapolated down to intersect the maximum-slope line (thick), 
giving the M2 value (large dot). Note: the quality of the 
Lorentzian fit is representative of some of the poorer fits. 
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FIG. 13: (a) Measured critical fields Hci,j and (b) locally 
saturated magnetizations Mj (full symbols) as a function of 
angle 9 measured from the [120] direction. Open symbols are 
reflections of the measured data across the 9 = 0° direction. 
Also shown are the calculated angular dependencies (dotted 
lines). 
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FIG. 14: Schematic representation of the three Ising-hke 
systems model (with three distinct R in the unit cell, sitting 
in unique orthorhombic sites): solid arrows- 'up', and dotted 
arrows- 'down' positions of the magnetic moments along the 
easy axes. 
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FIG. 15: Calculated magnetizations as a function of 6, based 

on the three eoplanar Ising model: open circles - S = 1 and 
crosses- S = 2 (see text for details). 
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FIG. 16: Polar plot ol the critical fields Hcij , with one possi- 
ble moment configuration shown for each observed metamag- 
netic state; open symbols represent reflections of the measured 
data -full symbols- across the 9 = 0° direction (see text). 
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FIG. 17: Anisotropic inverse susceptibility of TmAgGe (sym- 
bols) and the calculated polycrystallinc average (line); inset: 
low-temperature anisotropic susceptibilities. 
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FIG. 18: Anisotropic field-dependent magnetization, with 
the soUd hnes being the hnear fits of the high-field plateaus, 
extrapolated down to H = (see text). 




FIG. 19: M{e) of TmAgGc (full symbols) at T = 2 K and H 
= 70 kG, H_Lc. Solid line: calculated magnetization Mmax * 
cos{0 — n * 60°), n-integer. Note: = O" is defined at the 
[110] direction. 
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FIG. 20: (a) M(T) data for H = 1, 5 kG and 10 - 70 kG 
(AH = 10 kG) for H||[120], with enlarged high-field data in 
the inset; (b) M*T/H derivatives for H = 50, 60 and 70 kG, 

together with the Lorcntzian fits of the maxima (sohd linos), 
illustrating how the vertical line in Fig. 22c was determined 
for high applied fields. 
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FIG. 21: M(H) isotherms for T = 1.85, 2 - 4 K (AT = 0.25 
K), 4.2 and 4.6 K for (a) H||[120] and (b) H||[110]; insets 
show the enlarged data around the metamagnetic transitions. 
(Arrows indicate direction of increasing T.) 
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FIG. 22: H-T phase diagrams for TmAgGe, with (a) H||[110] 
and (c) H||[120], as determined from magnetization data in 
Fig.20-21; (b) intermediate-position phase diagram (see text); 
the error bar shown in inset of (b) represents the range of 
values for the high-field line, as determined from Fig.23. 
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FIG. 23: M(H) for T = 1.85 K (symbols), and the correspond- 
ing derivative (line). The inset shows an enlargement around 
the high field transition, to exemplify the two criteria used 
for determining this critical field (see text). Note: the error 
bar shown in the inset (determined from the position of the 
peaks in the derivative) gives a caliper of the uncertainty in 
determining this critical field value. 
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FIG. 24: M(H) isotherms (T = 2.0 K) for (a) O" < 61 < 10°, 
— 2" (enlarged Mi state shown in the inset) and (b) 
10° < 6 < 30", Ae = 2"; inset in (b) shows the extrap- 
olation of the two higher metamagnetic states down to lower 
fields, such that the intersection with the maximum-slope line 
gives the magnetization values M3 and Mcl-spm- solid dots 
(see text). Arrows indicate increasing 6. 
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FIG. 25: (a) Measured critical fields Hdj and (b) locally 
saturated magnetizations Mj (full symbols), as a function of 
angle measured from the easy axis. Open symbols are reflec- 
tions across the 6 = 0" direction (see text). Also shown are 
the calculated angular dependencies of Hdj and Mj (dotted 
lines). The error bars shown in the low part of (b) give the 
range of values that we can infer for Mi from the data shown 
in Fig.24. 
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FIG. 26: Polax plot of the critical fields He, with one of 
the possible moment configurations shown for each observed 
metamagnetic state. 



39 




FIG. 27: M(H) isotherms at T = 20 K for (a) 0° < 61 < 12° 

and (b) 12° < 6* < 30° {M = 1°); inset: enlarged M4 
state. Arrows indicate the direction of increasing 9. 
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FIG. 28: (a) Measured critical fields Hci,j and (b) locally 
saturated magnetizations Mj (full symbols), as a function of 
angle 8, for T — 20 K. Open symbols are reflections across 
the 6 = O'^ direction. Also shown are the calculated angular 
dependencies of Hdj and Mj (dotted lines). 




FIG. 29: Polar plot of the critical fields Hdj, with one possi- 
ble moment configuration shown for each observed metamag- 
nctic state (except for Mi, where the moment configuration 
is uncertain- see text); open symbols represent reflections of 
the measured data-full symbols- across the = 0" direction. 
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FIG. 30: Plausible models for extremely anisotropic, planar 
compounds (see text) with (a-b) tetragonal and (c-e) hexago- 
nal unit cells. The numbers are used to identify the different 
magnetic moments in the unit cell, (a) Four position clock 
model describing tetragonal systems with one R in tetragonal 
point symmetry; (b) Double coplanar Ising-like model: for or- 
thorhombic point symmetry in tetragonal unit cell, two mag- 
netic moments would be necessary (Ising-like systems, 90*^ 
away from each other in the basal plane), (c) Six position 
clock model: one magnetic moment with six possible orien- 
tations (arrows along six high symmetry orientations in the 
basal plane); (d) Triple coplanar Ising-like model: three R 
ions in unique orthorhombic point symmetry arc needed in a 
hexagonal compound (three Ising-like systems, 60" away from 
each other in the basal plane) ; (c) A double coplanar three po- 
sition clock model can describe hexagonal systems with two 
magnetic moments in unique trigonal point symmetry posi- 
tion, with three possible orientations (120° away from each 
other in the basal plane) for each. In all cases, the corre- 
sponding CL-SPM states (described in text) arc represented 
by full arrows (applied field is assumed to be vertically up. 



